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Abstract 

Design-oriented boundary conditions for subharmonic oscillations are of great interest recently. 
Based on a subharmonic oscillation boundary condition reported in a PhD thesis more than a 
decade ago, extended new boundary conditions are derived in closed forms for general switching 
DC-DC converters. Sampled-data and harmonic balance analyses are applied and generate equiv- 
alent results. It is shown that the equivalent series resistance causes the boundary conditions for 
voltage/current mode control to have similar forms. Some recently reported boundary conditions 
become special cases in view of the general boundary conditions derived. New Nyquist-like design- 
oriented plots are proposed to predict or prevent the occurrence of the subharmonic oscillation. 
The relation between the crossover frequency and the subharmonic oscillation is also analyzed. 
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1 Introduction 



Most power electronics textbooks [I] for example] devote two separate chapters on voltage mode 
control (VMC) and current mode control (CMC). For CMC, a separate slope-based design is used 
to avoid subharmonic oscillation. In switching DC-DC converters, three typical local instabilities 
[2] are subharmonic oscillation (period-doubling bifurcation, fast-scale instability), saddle-node 
bifurcation [3] (jump or multiple-solution instability), and Neimark-S acker bifurcation (slow-scale 
instability). Generally, most averaged models can predict the saddle-node bifurcation and the 
Neimark-S acker bifurcation [3], but not the subharmonic oscillation unless the sampling effects are 
considered as in [H [5J [B] . 

The subharmonic oscillation is known to occur when a sampled-data pole crosses -1 in the 
complex plane [2]. Some numerical boundary conditions of subharmonic oscillation based on sim- 
ulations have been reported for example]. Note that, here, the boundary condition means the 
critical condition in the converter parameter space, not about the well known critical eigenvalue at 
-1 in the complex plane. The boundary conditions define the subharmonic oscillation boundaries in 
the parameter space to separate stable and unstable regions. The numerical boundary conditions 
have limited usages for the converter design. Boundary conditions in closed-forms would greatly 
facilitate the converter design, because the quantitative effect of each relevant converter parameter 
can be clearly seen. Therefore, design-oriented boundary conditions for subharmonic oscillation are 
of great interest recently. For example, in [HH], the boundary condition for CMC with the voltage 
loop open is derived. In [9], a boundary condition in terms of ripple amplitude for a buck converter 
with a proportional-integral (PI) controller is derived. The PI controller used has only one pole 
and one zero. These boundary conditions do not consider the effect of equivalent series resistance 
(ESR). 

In [10], a unified VMC/CMC block diagram model is proposed, and a general closed- form 
boundary condition for the subharmonic oscillation is derived based on sampled-data analysis. 
Here, the term "unified" is used because the model is applicable to both VMC and CMC, and 
the term "general" is used because the boundary condition is applicable to all types of switching 
DC-DC converters of any system dimension under various control schemes. Although the general 
boundary condition was published in a PhD thesis []JJ] more than a decade ago, it was not published 
elsewhere and was not applied to various converters. It is reported here because some of recent 
results [U [TTJ associated with the subharmonic oscillation can be explained in terms of the general 
boundary condition. 

In this paper, based on the general boundary condition, extended design-oriented boundary 
conditions are derived. The boundary conditions in [JJ [H] for CMC with the voltage loop open and 
the boundary conditions in [9] for a buck converter with a PI controller become special cases in 
view of the general boundary condition. Additionally, the general boundary condition is extended 
for various converters under various control schemes. For example, it is extended to CMC with the 
voltage loop closed, average current mode control (ACMC) with a type II controller, and VMC with 
a type III controller. The type II or type III controller, popular in power electronics industry, has an 
integrator to improve the steady-state regulation and has more poles and zeros than a PI controller. 
The derived boundary conditions show the effects of various converter parameters, such as ESR, 
loading resistance, feedback gain, and pole/zero locations. It will be shown that the ESR causes the 
boundary conditions for VMC and CMC to have similar forms. For example, it is well known that 
the boundary condition for CMC has a term D — 1/2, where D is the duty cycle. It is also known 
[9] that the boundary condition for VMC has a term 1 — 2D + 2D 2 . It will be shown that, with 
ESR, the boundary condition for VMC or CMC has both terms of D - 1/2 and 1 - 2D + 2D 2 . The 
derived boundary conditions in this paper are expressed in terms of signal slopes, instead of signal 
amplitude as in 0. The slope-based boundary conditions agree coherently with the traditional 
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slope-based analysis for CMC. The boundary conditions can be also expressed in terms of other 
parameters, such as signal instantaneous slope, the source voltage or the compensating ramp slope, 
for example. It has been reported in |12| that the subharmonic oscillation is unrelated to the ripple 
amplitude for a converter under ACMC. It has also been reported in [13] that a converter can be 
stabilized by a compensating ramp with increased instantaneous slope but with similar amplitude. 

The analysis of subharmonic oscillation becomes simpler because of the derived closed-form 
boundary conditions. For a converter under a particular control scheme, once the converter is 
expressed in terms of the unified VMC/CMC block diagram model, the boundary condition for 
that particular control scheme can be readily obtained. For example, V 2 control |11] is of great 
interest recently because of its fast response. It will be shown that the boundary condition for V 2 
control can be easily obtained under the general modeling approach proposed in this paper. 

Although sampled-data analysis [IJJ |T5l |T6l |T7] has been known for many decades, the approach 
used here is different. First, traditional sampled-data analysis relies heavily on numerical analysis 
to obtain the steady-state solutions and the associated pole stabilities, which gives little insights 
for the converter design. In this paper, the emphasis is on the closed-form boundary conditions 
in the parameter space. The effect of each converter parameter can be clearly seen. Second, 
traditional sampled-data analysis generally starts from an approximate sampled-data model, where 
the accuracy would be lost at the beginning. In this paper, an exact sampled-data model is used 
to derive the boundary conditions which preserve the accuracy and are gradually simplified into 
various forms. 

Harmonic balance modeling \10\ [T8| fl~9| [20] is a complimentary approach to analyze subharmonic 
oscillations in switching converters |10[ 120]. The modeling approach proposed in |10[ 120] is further 
extended to analyze various control schemes and to derive similar boundary conditions as in the 
slope-based sampled-data analysis. This paper is originally planned in two parts (sampled-data 
slope-based analysis and harmonic balance analysis) because of its lengths. However, these two 
modeling approaches complement and corroborate each other. Both are presented here in a single 
paper for better cross reference. Note that this paper focuses on the sampled-data slope-based 
analysis. For some instances, harmonic balance analysis provides additional perspectives on the 
effects of some particular converter parameters. In this paper, the harmonic balance analysis is 
applied only to the buck converter. Similar analysis can be applied to other converters. Also note 
that this paper focuses only on the continuous conduction mode (CCM). 

The contributions along with the conclusion are stated at the end of the paper. The remainder 
of the paper is organized as follows. In Section [2] the operation of VMC and CMC is modeled 
in a single unified block diagram. In Section [3j steady-state analysis and small-signal analysis are 
presented. In Section HI general subharmonic oscillation boundary conditions for buck and boost 
converters are derived. In Sections [5] and [6] the sampled-data slope-based analysis and the harmonic 
balance analysis are applied respectively to analyze various control schemes. In Sections [7] and [HJ 
prediction of subharmonic oscillation based on the loop gain and crossover frequency are analyzed. 

2 Brief Review of Voltage/Current Mode Control Operation 

The operation of a general switching DC-DC converter in CCM under either VMC or CMC can 
be described exactly by a unified block diagram model |10[ [2T| [22] shown in Fig. [TJ A unique 
aspect about Fig. Q]is that the ramps in VMC and CMC are represented by the same hit). The 
operation of VMC and CMC is briefly reviewed here to make this paper self-contained. The control 
(reference) signal v r controls the output voltage v Q in VMC, and controls the peak inductor current 
i L in CMC. Denote the source voltage as v s . In the model, Ai,A 2 G K NxN , B 1: B 2 G R Nx2 , 
C,E 1 ,E 2 e R lxN , and D G R lx2 are constant matrices, where TV" is the system dimension. For 
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Figure 1: Unified VMC/CMC block diagram model for a switching converter in CCM 
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Figure 2: Waveforms for voltage mode control 



example, N = 5 for a buck converter with a type III compensator. Within a clock period T, the 
dynamics is switched between two stages, S\ and S 2 - Switching occurs when the ramp signal h(t) 
intersects with the compensator output y := Cx + Du € R. Denote the ramp amplitude as Vh, 
and denote the switching frequency as f s := 1/T and let uj s := 2irf s . 

Typical signal waveforms for VMC and CMC are shown in Fig. [2] and Fig. [3l respectively. 
In Fig. O the ramp has a positive slope, instead of a negative slope as commonly seen in most 
literatures, in order to be consistent with VMC. Other control schemes (average current mode 
control, for example) also fit the model of Fig. [TJ 



3 Brief Review of Steady- State and Small- Signal Analysis 

Steady-state analysis and small-signal analysis are briefly reviewed here to make this paper self- 
contained. The periodic solution x°(t) of the system in Fig. [1] corresponds to a fixed point x (0) in 
the sampled-data dynamics. A typical periodic solution x°(t) is shown in Fig. where x°(d~) = 
Aix°(d) + B\u and x°(d + ) = A 2 x°(d) + B 2 u denote the time derivative of x°(t) at t = d~ and d + , 
respectively. Let y°(t) = Cx°(t) + Du. In steady state, y (t) = Cx°(t). Let the steady-state duty 
cycle be D and d := DT. Confusion of notations for capacitance C and duty cycle D with the 
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Figure 3: Waveforms for current mode control 




Figure 4: A typical periodic solution x°(t) of a DC-DC converter in state space 



matrices C and D can be avoided from the context. 
In steady state, 



x°(d) = e Ald x°(0)+ / e Aia daB lU 

Jo 



x°(0) = e A *(T-d) x °( d ) + 



T-d 



e Ma doB 2 u 



(1) 

(2) 



From (P) and ([2]), one has 



x "(d) = (I-e Aia e 



A\d An(T- 



T-d 



e Ma daB 2 u+ / e AlCJ daB lU ) 



(3) 



Let B\ := [Sn,Si2], -B2 := [-E>2i,i?22] to expand the matrices into two columns. The buck 
converter generally has A\ = A2, B21 = 0at x 1i an d B12 = -622- For A\ and I — e AlT being 
invertible, (|3|) becomes 



x°(d) = (J - e AlT )^ 1 ^ 1 (e Ald - 1)^11^ - ^^12^ 



(4) 



Generally the controller may include an integrator (with a pole at zero), making Ai and / — e AlT 
non-invertible. In that case, the pole at zero can be replaced by a very small number 5, then A± 
and / — e AlT are invertible. Therefore, the invertibility of Ai or / — e AlT is not critical and can be 
resolved. This statement about invertibility of a matrix is not repeated later. 
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At the switching instant t = d, one has y°(d) = Cx°(d) + Du = h(d). Then, from 

h(d) + CA^ l B 12 v r - Du 
Vs ~ C{I - e A ^ T )- 1 A^ 1 {e A ^ d - I)B n 

The boost converter generally has B\ = B 2 , then 

rT-d r d 



(5) 



x°(d) = (I -e Ald e A2 ( T - d) )-He Ald [ e Ma da+f e Al ° da)B lU := X{d)B lU (6) 

Jo Jo 

Using a hat * to denote small perturbations (e.g., x n = x n — x°(0)), where x n is the sampled 
state at t = nT. From |10[ I21 |, [23] . the linearized sampled-data dynamics is 

x n+1 = <M n (7) 

where 

$ = e A 2 (T-d)(j _ (i°(d-) - x°(d + ))C Ald 
y*(d~)-h{d) 

Although the general methodology of sampled-data analysis has been known in the last three 
decades, the closed form expression of ([8]) was first published in [El [23], to the author's knowledge. 
The closed form of (|8|) greatly facilitates the derivation of the boundary conditions discussed next. 

4 Subharmonic Oscillation Boundary Conditions 

4.1 General Boundary Conditions 

The subharmonic oscillation occurs when one eigenvalue of $ is —1, and detfJ+c^] = 0. A necessary 
and sufficient boundary condition for occurrence of subharmonic oscillation is obtained |10[ p. 46], 

y°(d-) - C{e~ MiT - d) e- Ald + J)- 1 ^ ^) - x°(d + )) = h{d) (9) 
The proof is as follows. Suppose —1 is not an eigenvalue of e A 2(T- d ) e A id^ ^h en 

det[J - $] = det[I + e M{T-d) e A ld] det[/ _ (J + e A 2 (T- d ) e A ld) -i e A 2 (T- d ) ^) - i°(d + ) Ce A ld] 

Cx°(d~) — h(d) 

x°{d~)-x°(d+) 1 



y°{d-)-h{d) 

det[7 + <!>] = requires that the last term (inside the second square brackets) of the last equation 



det[J + e A ^ T ~ d ^e Md ][l - Ce Md {I + e MT-d) e A ld) -l & A 2 {T 

'equires that t 
equals to zero, which leads to Q. 

One can expand (|9|) in terms of x°(d), 

C(A lX °(d) + B lU ) - C(e- A2{T - d) e~ Ald + /) _1 ((Ai - A 2 )x°{d) + {B l - B 2 )u) = h(d) (10) 

Note that the condition ([9]) is valid for both VMC and CMC in this unified modeling approach, 
and it is applicable to general switching converters of any system dimension. Also note that in ([9]), 
the left side is related to the ripple slopes (y°(d~), x°(d~), and x°(d + )), and the right side is the 
ramp slope h(d). As in the popular slope-based boundary condition for the subharmonic oscillation 
in CMC, the condition ([9]) is also slope-based. 

The condition Q can be proved to be equivalent to 

y°(d + ) + C{e A2iT ~ d) e Ald + I)" 1 ^ ^") ~ x°{d+)) = h{d) (11) 



Using which one of (|9|) or (|lip depends on convenience. Since the proof for Q is given, the proof 
for ([TT]) is omitted to save space. 
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4.2 Buck Converter 



The buck converter generally has A\ = A2, B21 = Onxi, and Byi = i?22- Using (J4]), the boundary 
condition (I10p becomes 

C[(J - e AlT )- 1 (e Al<i -/) + (/ + e AlT )- 1 ] J Bn = ^ (12) 
or in terms of v s , which shows the critical value of v s , 

hid) , \ 

s C[(/-e^ T )- 1 (e Ald -/) + (/ + e^ T )- 1 ] J Bn { } 

Subharmonic oscillation is avoided if v s < V* (if the denominator of (fT3|) is positive as discussed 
later). The condition f)13j) is for the trailing edge modulation [23] . For the leading edge modulation, 
a similar condition has been reported in [10, p. 72]. 

As will be shown later for the VMC buck converter, the curve (|13p as a function of d = DT 
generally has a minimum at d = T (or equivalently D = 1) with a value 

V*\ = V*\ = C\A\ 

s lmm sU=T C[{I + e A ^ T )~ l - I]B n 1 ; 

Then, for v s < V*\ min , the subharmonic oscillation is avoided for all duty cycles. 

Based on the assumption that the switching frequency f s = 1/T is much larger than the 
absolute value of any eigenvalue of A±, matrix approximations such as e AlT ~ 7 + A\T + A\T 2 /2 
and (I + AiT)^ 1 ml — A\T can be applied. Then, the boundary condition (|12p leads to 

^.(^w.M (15) 

2 4 v., 



or in terms of V^* , 



V* » ^ 5 (16) 



Throughout the paper, all approximate boundary conditions are expressed with the approximation 
sign and all exact boundary conditions are expressed with the equality sign "=." 
Remarks: 

(a) The left side of (|15p is a weighted combination of CBn and CAiBnT. It will be shown 
that if the equivalent series resistance (ESR) R c = 0, CA\B\i dominates in VMC, whereas CB\\ 
dominates in CMC. For R c > 0, either VMC or CMC has both the terms CB n and CA\B X \T, 
indicating that the subharmonic oscillation conditions for VMC and CMC are closely related. This 
result has not been reported. 

(b) The boundary condition (|15p seems to have a pattern. It is a hypothesis that the exact 
boundary condition has the form 



C( V 6 n (D)A?T n )B n = ^ (17) 

^ ■* ?!. 

71=0 



where 5 (D) = (1 - 2D)/2, 61(D) = (-1 + 2D - 2D 2 )/4, etc., and d8 n+1 (D)/dD = 5 n (D). Since 
an exact condition as in (|12p has been obtained, another exact condition in series expression may 
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give additional insights but may be unnecessary. As a side note, 62(D) = (-D + 3D 2 — 2D 3 )/12 
and 5 3 (D) = (1 - 2D 2 + 4L> 3 - 2L> 4 )/48. Note that d6 2 (D)/dD / 61(D), but d5 3 (D)/dD = 5 2 (D). 
Further research on the series expression of (|17p is pursued and will be reported separately. 

(c) The condition (|12|) is an exact condition, whereas (I15|) is an approximate one. Based on 
simulations, if the real controller poles are smaller than w s /10, the approximate condition (|15p 
is close to the exact condition (|12p . and (|15p is generally adequate to predict the subharmonic 
oscillation. If the real controller poles are greater than u s /10, using the exact condition (j 12 j) gives 
more accurate results. 

(d) The boundary condition (|15p does not require a matrix inverse. 



4.3 Boost Converter 

The analysis for the boost converter is similar to that for the buck converter. Let A(d) := I + 

(Ai - (I + e~ A2 ^ T ~ d ' ) e~ Ald )~ l (Ai - A 2 ))X(d) to simplify the equation. Using ©, the boundary 
condition ([TU]) becomes 

CA{d)B lU = h{d) (18) 

or in terms of v s , 

Vs ~ CA(d)B u (19) 

Based on (|19p. it can be proved that, in a boost converter in CCM with proportional feedback VMC 
(thus the system dimension N = 2) and with practical converter parameters, the subharmonic os- 
cillation does not occur. With different modes or control schemes (with a higher system dimension), 
subharmonic oscillations may still occur in a boost converter. For example, subharmonic oscillation 
occurs in a boost converter in discontinuous conduction mode 1251. 



4.4 The "S plot": a Slope-Based Plot in the Real Domain 

Define an "S plot" as a function of D = d/T, 

S{D) := y°{d-)-C{e- A ^ T - d ^e- Ald + I)- 1 {x {d-)-±°{d + )) (20) 

f C[(J - e AlT )^ 1 (e Ald -/) + (/+ e^ 171 )- 1 ]^!!^ (for the buck converter, from (fT2|\ 
~~ \ CK(DT)Biu (for the boost converter, from JIB])) ^ ' 

Then, from Q, subharmonic oscillation occurs when 

S(D) = h(d) (22) 

The S plot facilitates the converter design to avoid the subharmonic oscillation. For example, if 
the ramp slope h(d) is large enough such that h(d) > S(D), then the subharmonic oscillation is 
avoided. 

Note that, for the buck converter, S(T) = C[(I + e A ^ T y x - I]e Ald B u v s w -(CB n /2 + 
CA]_B\\T/£)v s . It will be shown later that, for the VMC buck converter, the maximum of 
S(D) is generally S(T). To avoid the subharmonic oscillation for all duty cycle, a ramp slope 
h(d) > max S(D) is required. It may be over-compensating because such a ramp is applied to 
avoid the subharmonic oscillation for all duty cycle. To avoid the subharmonic oscillation for a 
particular duty cycle, only h(d) > S(D) is required. 
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4.5 Other Approximate Slope-Based Boundary Conditions 

In the analysis above, a (first) approximate slope-based boundary condition (fT5j) is presented. 
Other approximation approaches may be applied to derive similar slope-based boundary conditions. 
Generally these boundary conditions are close to the first approximate boundary condition (|15p . 

4.5.1 Second Approximate Boundary Condition 

With a large switching frequency, (I + e M{T-d) ^dyi _ j j 2 _ + a 2 {T - d))/4. Then the 
boundary condition ([TT]) becomes 

\{y\d-) + y\d + )) - j(A,d + A 2 (T - d))(x (<T) - x°(d+)) « (23) 

This condition is also expressed in terms of signal slopes. Then, a similar condition as in (|15p can 
be obtained. 

4.5.2 Third Approximate Boundary Condition under a Large Switching Frequency 

Although the switching frequency never reaches infinity, the results obtained under this condition 
link very well with the well known stability condition for CMC. Note that this section is presented 
here only to illustrate a special case for the boundary condition ([9]) . When the switching frequency 
is high, e Ald ~ e A ^( T ~ d ) ~ J. Then the boundary condition becomes 

±(f(d-) + y (d+))*h(d) (24) 

Like (fT2|) . (}2"4"]) is also a slope-based condition. In CMC, the (inductor current) slopes y°(d~) 
and y°(d + ) and the ramp slope h(d) are generally expressed in most textbooks [U p. 448] as —mi, 
rn,2, and m a , respectively. Then, (|24p corresponds exactly to the well known minimum ramp slope 
m a = (777,2 — mi)/2 required to stabilize the converter. 

As a side note, the condition ([9]) or (|24p occurs when $ has an eigenvalue —1. Similar to the 
proof for ©, one can prove that y°(d + ) = h(d) when $ has an eigenvalue 0, which generally causes 
a deadbeat effect. Under CMC, this condition is equivalent to the well known condition ni2 = m a 
to have a deadbeat effect pp. Note that y°(d + ) = h(d) is a general condition valid for any control 
scheme, not just CMC. 

5 Application of Sampled-Data Slope-Based Analysis 

Without loss of generality, the sampled-data slope-based analysis is applied to buck converters. 
Analysis of Boost converters can be applied similarly. The exact condition (|12p and the approx- 
imate condition (|15p are applied to analyze various VMC/CMC control schemes presented next. 
The compensator complexity increases from a simple proportional compensator to a type III com- 
pensator. The system dimension iV increases from two to five. 

5.1 Proportional Voltage Mode Control (PVMC) 

Without loss of generality, let the voltage loop have a proportional feedback gain k p . One has y = 
k p (v r — v ), as shown in Fig. For designation purpose, this control scheme is called proportional 
VMC (PVMC). 
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Figure 5: A buck converter under PVMC 



In the power stage, let the state be x = (iL,vc)', where %l is the inductor current and vq is 
the capacitor voltage. Let the load be R, the inductance be L, the capacitance be C, and the 
equivalent series resistance (ESR) be R c . Let p = R/(R + R c ). For R c = 0, p = 1. Then, 



Bi 
Ex 



P 



-Rc 
L 
1 
C 



[Bu,Bi2\ = 
E 2 =p[R c 



-l 

L 

-1 





[B 2 \ , B 



22J 








(25) 

(26) 
(27) 



Then, from (d§D, 



For i? c = 0, 



C = - 
For R, 



V* 



k p p[R c , 1] 

_ -kppR c 



L 

0, CBn 



D = [0, 



(1 



R 2 C C , 
L ' 



hp 
TC 



4V h LC 

pk p T* { MkC( D _ i) + p (i _ _ 2D + 2^) ' 



/jfcpT 2 y - \) + {I - 2D + 2D 2 ) 



) 



1 



1 - 2L» + 2L» 2 



(28) 
(29) 

(30) 



In [9], a similar condition (only for R c = 0) with a minor correction term is obtained. Compared 
with ()30p . (|28p has an additional term involving D — 1/2, which is related to CMC when the 
subharmonic oscillation occurs as discussed next. For R C C <C T, (|30|) is a good estimate of the 
critical source voltage. If the condition R C C <C T is not met, (j28[) will be more accurate than ()30p . 

Both (|28|) and (|30p are functions of D, and D is chosen to meet the required output voltage. 
For example, v Q ~ for the open-loop buck converter. Given a value of D, (|28p or (|30p sets 
the maximum value of t> s to avoid the subharmonic oscillation. In a closed loop converter, D is 



11 




0.2 0.4 0.6 o.e 



Duty cycle 

(a) 




(b) 



Figure 6: Plots of (|28|) (solid line) and (|31|) (dashed line), and the intersection shows V*. (a) 
R c = 0, (b) R c = 2 mO. 



determined by ([5]), or by the steady-state equation, DVh = k p (v r — v Q ) = k p (v r — v s D), rearranged 
as 

Vs = d~T p (31) 

For a practical converter, k p is large and v s ~ v r /D. When the subharmonic oscillation occurs, 
both the steady-state condition f)31f) and the subharmonic oscillation condition ()28f) (or (|30p ) need 
to be met. Therefore, the intersection of (I28p (or (|30f> ) with (I3ip is the critical source voltage in a 
closed loop converter. 

Another way without plotting the two curves (|28p and (|3ip to determine the critical value for 
the subharmonic oscillation is as follows. One way is to subtract (I3ip by (I28p . Another way is to 
divide (f3T|) by ([28]) . which leads to an "M plot" as discussed in Seed 

Example 1. Consider a PVMC buck converter with k p = 80. The converter parameters are 
V h = l,v r = 4,f s = l MHz, L = 1 juH, C = 100 fj,F, and R = 2 ft. 

To see the effects of R c , the curves of (|28p and (|3ip for different values of R c are also shown in 
Fig.EJ The intersection of (|28p and (|3ip determines V,*. The converter is stable for those operating 
range of D such that the curve (f3T|) is below the curve (f28"p . It can be shown that the curve ([28]) 
is actually almost identical to the curve (|13|) . omitted for brevity. 



For R c = 0, the curves of (|28p and (|3ip are shown in Fig. 6(a), intersecting at (D,v s ) = 
(0.41,9.7). The stable operating range of D is [0.41,1]. Simulation (Fig. [7]) with v s = 10 > V* 
shows the subharmonic oscillation as predicted. 

For R c = 2 mO, the curves of (I28p and (I3ip are shown in Fig. |6(b)[ intersecting at two points 
(D,v s ) = (0.34,11.85) and (0.89,4.48). The curve ([28D is shifted to the upper-left (because the 
term (4R C C/T)(D — 1/2) in (|28p becomes significant), and V* is increased to 11.85. Simulation 
(Fig. [8]) with v s = 12.5 > V* also shows the subharmonic oscillation as predicted. The stable 
operating range of D is [0.34, 0.89]. For v s = 4.44 (and D = 0.89), another subharmonic oscillation 
occurs, along with a border-collision bifurcation, because y(t) is out of bounds of h(t). The signal 
plot is omitted for brevity. □ 
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Figure 7: Plots of y(t) and ix(t) (solid lines), and h(t) (dashed line), R c = 




Figure 8: Plots of y(t) and ii{t) (solid lines), and h(t) (dashed line), R c = 2 mO 



5.2 l/ 2 Control: Similar to PVMC 

The constant-frequency peak voltage regulator (CF-PVR, as shown in Fig. [9]), a type of V 2 control 
[llj . is proved below to be a special case of PVMC in terms of the boundary condition. In CF-PVR, 
the output voltage is sensed (through a voltage divider with a gain of k p ), added with a stabilization 
ramp h(t), and compared with a reference signal v r to determine the duty cycle. In terms of Fig. [H 
y = v r — k p v , and the matrices C = —k p p[R c , 1] (same as for PVMC) and D = [0, 1]. 

Since, in terms of Fig. [H the model is almost the same as PVMC, the boundary condition (|29p 
for the PVMC buck converter is also applicable to CF-PVR. Assume R c R, then p ~ 1. From 
(|29p . the subharmonic oscillation is avoided if 

. WhLC 1 

Vs< k p T* [ ^Ml(D- \) + {\-2D + 2D*y 1 ' 

With the fact that v Q = Dv s for the buck converter, one can rearrange (|32p as 

V h L 2D -I T A — 2D D . 

>^T^ + lT^(—^ + ^) (33) 



k pVo R c T 2D R c C y AD 2 
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Figure 9: A buck converter under CF-PVR 

which agrees with [HI Eq. 7] (by setting k p = 1) and shows the required ramp amplitude Vh to 
stabilize the subharmonic oscillation. Without the ramp (Vh = 0), rearranging (I32p . subharmonic 
oscillation is avoided if 



2 ^ 1-2D 



or, equivalently, 



R C C 1 D 2 1 

— >2 + T=W* D< 2 (35) 



also agreed with [111 Eq. 5]. 
Remarks: 

(a) The conditions (|33p is applicable to &ot/i PVMC and CF-PVR, and (|34p or (|35p is a special 
case (T4 = 0) of ([33|). 

(b) .D < 1/2 is explicitly required in (|35p . whereas D < 1/2 is implicitly required in ()34l) . 

(c) The subharmonic oscillation for CF-PVR does not necessarily occur at D = 1/2 as claimed 
in [26]. Based on (|35p . if i? c C = 5T/8, for example, the subharmonic oscillation occur at D = 1/4. 
A single pole is shown in [261 P- 352] to be unstable when D > 1/2, it cannot be a subharmonic 
oscillation because a pair of unstable complex poles is required. 

5.3 CMC With the Voltage Loop Open 

In CMC with the voltage loop open, as shown in Fig. [TUl a switching occurs when i c — ii(t) = h(t), 
where i c (denoted as v r in Fig. [1]) is a peak inductor current control signal. One has y = —ii + i c 
and 

C =[-1,0] D = [0,1] 

CB U = =jk CA 1 B 11 = 

For R c = 0, CB U = =± CA x B xl = 

Note that in CMC, CB\\ dominates (much greater than CA\B\\F in absolute value). Let the ramp 
slope be m a = h(d), and let the critical (minimum) ramp slope required to stabilize the converter 
be m*. 

Case I: R r = 0. 
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Figure 10: A buck converter under current-mode control with the voltage loop open 

For R c = 0, from (|15p . subharmonic oscillation is avoided [8] if 

m a> m* a ^ V -±(D- 1 -) (36) 

Without the compensating ramp (m a = 0), the converter is unstable for D > 0.5 as the condition 
(|36p is violated (because the left side of (j36|) is zero whereas the right side of (j36|) is positive). The 
condition (136p defines the minimum of m a and L, and the maximum of v s and D to avoid the 
subharmonic oscillation. For example, the maximum of v s to avoid the subharmonic oscillation is 
m a L/(D — 1/2), which shows that the stable operation range of the source voltage v s is linearly 
proportional to the ramp slope. Adding the compensating ramp not only improves the stability 
but also enlarges the operating range of the source voltage (line regulation). 

Case II: R c > 0. 

For CMC with R c > 0, the boundary condition (|15p becomes 



m a ™ T (D--- — ( )) (37) 

Compared with (|36p . (|37p has an additional term 1 — 2D + 2D 2 related to VMC as discussed above. 
Without the compensating ramp (m a = 0), ([37]) is a quadratic equation of D. Ignoring the second 
solution greater than one, the critical value of D when the subharmonic oscillation occurs is 



D * „ i ± L . L R2 c T2 

2 R C T R C T V 4L 2 

W \ + E ^f {i ° xRcT ^ L) (38) 

The subharmonic oscillation is generally believed to occur when D = 1/2 without the ramp. 
However, from ()38|) . the ESR extends the operating range of the duty cycle beyond 1/2. Generally, 
pR c T/8L is small, and the instability at D = 1/2 is still a good criterion. 

The boundary conditions can be also rearranged in terms of v s . From (|36|) and (|37|) . for R c = 
and R c > 0, respectively, 

V? ~ ^(— ^0 (39) 
« YhLt ) (40) 

^2 L \ 4 / 
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Figure 11: A buck converter under current-mode control with the voltage loop closed 
Remarks: 

(a) Care should be taken about the direction of instability. In an equality, when the both sides 
are divided by a negative number, the inequality sign is reversed. Assume that the ramp is very 
small, for D > 1/2, the converter is stable if v s < V*. For D < 1/2, the converter is stable if 
v s > V*. 

(b) Comparing these two equations with (|28|) and (|3U|) for VMC, one can see that the boundary 
conditions for VMC and CMC have very similar forms. Both have the terms D — 1/2 and 1 — 2D + 
2D 2 for R c > 0. 

5.4 CMC With the Voltage Loop Closed 

Without loss of generality, let the voltage loop have proportional feedback k p (v r — v Q ), as shown in 
Fig. CCD One has y = k p (v r — v ) — ii and 

C = -[l + k pP R c , k pP ] D = [0,k p ] 

CB n = =±{k pP R c + 1) CAxBn = 4^(1 - ^) + ^ 

For R c = 0, CBu = =± CAiB u = ^ 

The boundary condition (fT5j) rearranged in terms of the critical feedback gain k p is 

m a L , pR c T ( 1-2D+2D 2 \ _n , 1 

k * ~ «■ + L [ * ] + 2 C411 

For a special case i? c = 0, one has 

m a L _ j-j , 1 

P T ( 1~2D+2D 2 \ ^ > 

C'\ 4 ^ 

agreed with |27t Eq. 30]. A larger T/C leads to a smaller fc*. Also, as D is closer to 1/2, 
(1 - 2D) /(I - 2D + 2D 2 ) and hence fe* become smaller. 

Subharmonic oscillation is avoided if k p < k*. The gain margin to avoid the subharmonic 
oscillation is 201og(£;p). The gain margin based on traditional averaged analysis generally helps to 
avoid the Neimark-S acker bifurcation, but not the subharmonic oscillation, as shown in the next 
example. 
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The traditional CMC design is as follows. Choose a ramp slope to ensure the current loop 
is stable. Based on an averaged model, one derives the control-to-output transfer function and 
designs a voltage loop controller to have enough phase and gain margins. It is possible that even 
with current loop designed to be stabilized by the compensating ramp, the subharmonic oscillation 
still occurs due to the voltage loop. 

Example 2. (Slope-based analysis can predict the critical value of feedback gain k* to avoid the 
subharmonic oscillation, whereas the averaged models fail to predict it.) Consider a CMC buck 
converter with proportional control (compensator transfer function G c (s) = k p ) in the voltage 
loop. The converter parameters are f s = 300 kHz, L = 900 nH, C = 990 /jF, R = 0.4 £1, and 
v r = 3.34 V. The source voltage is 5.5 V, and the duty cycle D is around 0.6. Since D > 0.5, 
the converter is unstable if a compensating ramp is not added. Adding a compensating ramp 
m a = 7712/2 = v s D/2L = 1.8333 x 10 6 V/s is supposed to stabilize the converter for all duty cycles, 
where ni2 is the inductor current slope during the off stage. One has Vh = = 6.11. However, 
with the voltage loop closed, the subharmonic oscillation may still occur even though the current 
loop is supposed to be stabilized by the ramp. 

The ESR R c contributes to the output voltage ripple, and hence affects the current loop. Two 
cases are simulated: R c = 5 mil and R c = 0. Simulations show that with R c = 5 mfi, the critical 
gain k* is 237, and with R c = 0, k* = 452, which are also confirmed by the exact sampled-data 
analysis [21] . 

Take R c = 5 vaD, and k p = 237, for example. The signal waveforms indicating (weak) subhar- 
monic oscillation are shown in Fig. [12j The steady states are 2T-periodic orbits. The unstable 
T-periodic orbit y°(t) = Cx°(t) + Du = k p (v r — v Q (t)) — Il{P) (with il inverted) still exists and 
is shown in Fig. [T3l along with h(t). If the T-periodic orbit x°(t) is perturbed, it will go to the 
stable 2T-periodic orbit shown in Fig. [T2J The designed slopes of (on-time inductor current slope) 
mi = (1 — D)v s /L and (off-time inductor current slope) = Dv s /L are 2.48 x 10 6 and 3.63 x 10 6 , 
respectively. However, ESR R c contributes to the output voltage ripple, and the actual slopes of 
y = k p (v r — v ) — iL shown in Fig. [T3]are 5.39 x 10 6 and 7.88 x 10 6 , respectively. Even with increased 
ripples, (y°(d + ) + m a )/(y°(d~) + m a ) = —0.84 > —1, and the current loop should be still stable 
according to the averaged models. However, the subharmonic oscillation still occurs as shown in 
Fig. m 

The critical gain k* can be predicted by (|4"Tj). From (|4"T]) . k* = 223 for R c = 5 mf2, and k* = 468 
for R c = 0, agreed closely with the simulation results and the exact sampled-data analysis. The 
small discrepancy is due to the error of the duty cycle. With the voltage loop closed, the actual 
value of the duty cycle from the simulation for R c = 5 mfi is 0.5941, not 0.6 as supposed. With 
the actual value of the duty cycle of 0.5941, from one has k* = 237, agreed exactly with the 
simulation results. 

Two averaged models are used for comparison. One is from [3], and the other is from [TJ p. 
470]. The control-to-output (i c -to-v Q ) frequency responses are, respectively, also shown in Fig. [Til 

4441321980(s + 2.02 x 10 5 ) 
(s + 3273)(s 2 + 5.92 x 10 5 s + 8.883 x 10 11 ) 

5000(s + 2.02 x 10 5 ) 
(s + 3276)(s + 9.993 x 10 5 ) 

Both show infinite gain margins, and do not accurately predict the critical gain of 237 to avoid the 
subharmonic oscillation. Both models have a similar ESR zero and a similar low frequency pole, 
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(43) 



whereas the model of [1] has a pair of complex poles associated with half the switching frequency 
as the resonance frequency instead of a single pole as in [I]. 

With k p = 237, the loop gains (i c -to-v k p ) of the two averaged models are shown in Fig. [15j 
The averaged model of [4] shows that the phase margin is 36.5 degrees. However, the subharmonic 
oscillation still occurs. The crossover frequency is w c = 1.28 x 10 6 rad/s, and w c /uj s = 0.68. It will 
be shown in Section [5] that a large crossover frequency likely leads to the subharmonic oscillation. 

The averaged model of [1] even shows a phase margin of 129 degrees. However, the subharmonic 
oscillation still occurs. The crossover frequency is w c = 7.19 x 10 5 rad/s, and w c /oj s = 0.3814 is 
also large as will be discussed in Section [H □ 




0.2 0.4 0.6 0.8 1 1.2 

Time (Second) 



Figure 12: Signal waveforms showing (weak) subharmonic oscillation, R c = 5 mQ and k p = 237 
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Figure 13: Unstable orbit y°(t) (solid line) and h(t) (dashed line), R c = 5 m£l and k p = 237 

5.5 Enhanced V 2 Control: Similar to CMC With the Voltage Loop Closed 

Compared with CF-PVR, the enhanced V 2 control [26] has an additional loop from the inductor 
current through a sensing resistor Ri shown in Fig. [TBI One has y = v r — v — R^l- Similar to 
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(a) (b) 
Figure 14: Bode plots for (a) averaged model of [3], (b) averaged model of [JJ 



Bode Diagram 

Gm = Inf dB (at Inf rad/sec) , Pm = 36.5 deg fal 1 .28e+0D6 rad/sec) 



Bode Diagram 
Gm = Inf , Pm = 1 29 deg (at 7.1 9e+005 rad/sec) 





Frequency (rad/sec) 



Frequency (rad/sec) 



(a) (b) 
Figure 15: Loop gains for (a) averaged model of [3J and (b) averaged model of pQ 



CMC with the voltage loop closed and k p = 1, one has 



C = -[pR c + Ri,p] 
CB n = =£(pRc + Ri) 



£ = [0,1] 
CAxBu = 



=£1(1 

LC V 1 



pL ' 



Consider a special set of converter parameters (with large capacitance C and small ESR R c ) \26\ 
p. 351], for example. If L/R c > (R c + Ri)C > T, then CAiBn > CA\BuT and the subharmonic 
oscillation occurs (without the ramp, h(t) = 0) around D = 1/2, agreed with \2Q\ p. 351]. Similarly 
to the V 2 control, the subharmonic oscillation does not occur exactly at D = 1/2. For other set of 
converter parameters, the analysis is similar to CMC with the voltage loop closed and is omitted 
here to save space. 
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Figure 16: A buck converter under enhanced V 2 control 




Figure 17: System diagram of an ACMC DC-DC converter 



5.6 Average Current Mode Control (ACMC) with a Type II Compensator 



An ACMC DC-DC converter is shown in Fig. El The operation of ACMC is as follows [25]: The 
inductor current is sensed by a resistor R s and compared with a voltage reference v r from a 
voltage loop (not shown). The difference is amplified by a current-loop compensator, generally a 
type II compensator [29J P- 256], 



G c (s) :-- 



H n (s) K c (l + ^) 
H d (s) ( s + 6)(l + - 



(45) 



uj p , where 5 is a small number close to zero. For 5 = 0, 



The compensator has poles at —5 and 
the compensator has an integrator with a pole at zero. The compensator output y(t) is compared 
with the ramp h(t) to determine the duty cycle. Let the state be x = (il, vc, v e i, v e 2)', where v e \ 
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and v e 2 are the states of the current-loop compensator. Then, 
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Since the integrator is modeled as l/(s + 5), both A\ and A 2 are invertible. 
Case I: ui p > u s /lQ. 

In this case, the compensator pole oj p is close to oj s , and those higher order terms (CA 2 BnT 2 , for 
example) in (|17p may be larger than CBn and CA\BnT, and one needs to use the exact condition 



AiT\-l 



T3|) . Generally, the compensator has a pole at zero in the complex plane, making (/ — e 1 ) 



singular in (|13p . It can be easily resolved by using a very small number for the zero pole as discussed 
previously. 



Case II: ui p < u; s /10. 

Based on simulations as shown in the next example, for ui p < u s /10, the boundary condition (|16p 
is still a good approximation. Here, CB\\ = and CA\Bn = —K c R s u> p /uj z L, then the boundary 
condition (|16p becomes 

4W 1 

s T 2 K c R s Lo p { l-2D + 2D 2) { } 



Example 3. Consider a buck converter under ACMC [30\ p. 114]. The power stage parameters 
are v s = 14 V, v a = 5 V, v r = 0.5, Vh = 1, /, = 50 kHz, L = 46.1 pR, C = 380 pF with ESR 
R c = 0.02 f2, and R = 1 £1. The inductor current sensing resistance is R s =0.1 The compensator 
has a zero at ui z = 5652.9 rad/s, two poles at and u p , and a gain if c = 75506. 

The compensator pole oj p is varied from 0.14w s to 0.81w s . An unstable window of ui p between 
0.18a; s and 0.49^ was found and reported in [12]. When ui p is inside the window, the subharmonic 
oscillation occurs. Take ui p = 0.49o; s , for example. The converter is unstable. The signal waveform 
indicating the subharmonic oscillation is shown in Fig. [TSl The unstable T-periodic orbit is shown 
in Fig. fl9l When the T-periodic orbit is perturbed, it will lead to the 2T-periodic orbit shown in 
Fig. [HJ The unstable window of ui p can be predicted exactly by (|13p . shown in Fig. [20j 

In the above analysis, v s is fixed at 14 V, and ui p is varied. Next, let ui p be fixed at w s /10, and 
v s is varied to determine V* . Since uj p is smaller, the approximate condition (|46p is close to the 
exact condition (|13p as shown in Fig. [2TI indicating V^* = 19. From (|46|) . a smaller ui v results in a 
larger V^* as expected. 

In ACMC, the subharmonic oscillation is unrelated to the ripple amplitude |12j . Take ui p = 
0.81u; s , for example. The converter is stable. The stable T-periodic orbit is shown in Fig. [221 In 
Fig- E21 the signal y°(t) has a large signal amplitude but is still stable, whereas in Fig. \W[ y°(t) 
has a smaller signal amplitude but is unstable. □ 
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5.7 Average Current Mode Control (ACMC) with a PI Compensator 

Instead of a type II compensator, consider an ACMC buck converter with a PI compensator, 
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(47) 
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Figure 20: Plot of (|13p . intersecting with v s = 14 at oj„/u s = 0.18 and 0.49, indicating the unstable 
window of u p 
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Figure 21: Approximate condition (|46p (dash-dotted line) is close to the exact condition (|13p (solid 
line), V,* = 19 



Since the integrator is modeled as l/(s + 5), both and A2 are invertible. One has (with 5 ~ 0) 

CB U = -*^, CA lBll = ^(^f-l) (48) 

Compared with the type II compensator where a high frequency pole is included, the PI compen- 
sator does not have a high frequency pole. The boundary condition (|15p . expected to be accurate 
without the high frequency pole, becomes 

2 4 L v s K s K c 

which agrees with the boundary condition derived in [31|, Eq. 4.6] (based on the describing function 
method where the effect of R c is ignored and pR c <C uj z L is assumed) The boundary condition is 
rearranged as 

1 A - 2D + 2D 2 , h(d)Lu) z , , 
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Figure 22: Stable T-periodic y°(t) with a larger ripple, oj p = 0.81w s 



Power Stage 




Figure 23: A DC-DC converter with a type III compensator G c (s) = H n (s) / Hd(s) 

For T{ui z — pR c / L) <C 1 and from (l4*9j) . the ramp slope required to avoid the subharmonic oscillation 
is 

w>^0 < 51 > 



5.8 VMC with a Type III Compensator 

A type III compensator \29\ p. 261] has three poles, two zeros, and a gain K c , with a transfer 
function 

H n (s) _ K c (l + ±)(l + ±) 



G c (s) 



H d (s) ( a + (J )(i + i)(l + i) 



(52) 



where 5 is a small number close to zero. For 5 = 0, the compensator has an integrator with a pole 
at zero. 

For a buck converter with a type III compensator as shown in Fig. [23l let the state be x = 
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(*L? vc, v e ±, v e 2, v e 3)', where v e \, v e 2, and v e s are the states of the type III compensator. Then, 
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Since the integrator is modeled as l/(s + <5), both Ai and A 2 are invertible. 
Case I: p\ > w s /10. 

A typical guideline [291 P- 412] popular in industry to set the parameters of the compensator is as 
follows. Set one pole at 6 ~ (as an integrator), and set pi = oj s /2 and P2 = 1/R C C. Set the gain 
K c to adjust the phase margin and the crossover frequency. Let z\ = K z j\fLC and z 2 = l/^/LC, 
where k z is a zero scale factor to have additional flexibility to adjust the phase margin and the 
crossover frequency. The zero scale factor k z used in industry typically varies between 0.1 and 1.2. 
As will be shown later, a smaller value of k z may lead to the subharmonic oscillation. Taking into 
account the above guidelines, the compensator has a transfer function 



The next example shows that, with the compensator (|53p. the subharmonic oscillation still 
occurs even with a phase margin of 38.9 degrees. As in ACMC, one compensator pole is close to 
cj 8 , and those higher order terms (CAfBuT 2 , for example) in (|17p may be larger than CBn and 
CAiBuT, and one needs to use the exact condition (fT3j) , 

Example 4. (With phase margin of 38.9 degrees, the subharmonic oscillation still occurs.) 
Consider a buck converter with the type III compensator ()53|) . Exactly the same parameters as in 
[32] are used: f s = l/T = 300 kHz, L = 900 nH, C = 990 /uF, R = 0.4 SI, R c = 5 mfi, d s = 5V, 
v r = 3.3 V, V h = 1.5 V, K c = 7.78 x 10 4 , z x = 1/2VLC = 1.675 x 10 4 , z 2 = 1/VIC = 3.35 x 10 4 , 
Pl = uj s /2 = 9.425 x 10 5 , and p 2 = 1/R C C = 2.02 x 10 5 . 

Simulation (Fig. I24j) shows that the subharmonic oscillation occurs when v s = 16 V (D ~ 0.206). 
This is also confirmed by the exact sampled-data analysis |21j with a sampled-data pole at -1 when 
the subharmonic oscillation occurs. The loop gain frequency response (Fig. 1251) shows that even 
with a phase margin of 38.9 degrees, the subharmonic oscillation still occurs. 

The subharmonic oscillation can be accurately predicted by the slope-based analysis. If D 
is designed or known as 0.2, for example, then V^* = 15.6 based on (|13p . which is close to the 
simulation results and the exact sampled-data analysis. Accuracy can be improved by using (|13p 
and ()3ip to determine the right duty cycle, shown in Fig. [26"1 which shows exactly that V* = 16 
and D = 0.206. □ 




(53) 



(s + 5)(l + ^)(l + R c Cs) 
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Figure 24: Signal waveforms showing the subharmonic oscillation 
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Figure 25: Loop gain frequency response showing phase margin of 38.9 degrees, but the subharmonic 
oscillation still occurs. 



Example 5. (Unstable window of p\, unrelated to the ripple size of y°(t).) Consider again 
Example [H Instead of v s = 5, let v s = 16, the subharmonic oscillation occurs. The stable 2T- 
periodic orbit y°(t) is similar to Fig. 1241 omitted to save space. The unstable T-periodic orbit y°(t) 
is shown in Fig. [271 

Here, p\ = 0.5w s . In the following, p\ is varied from 0.1w s to 0.6w s to see the effect of this 
pole location on stability. Similar to Example [3l the value of p\ adjusts the ripple size of y°(t). 
A larger p\ leads to a larger ripple of y°(t). In [9], it is hypothesized that the ripple size of 
y°(t) is related to subharmonic oscillation. The following simulation shows that the ripple size of 
y°(t) is unrelated to subharmonic oscillation. The S plot (Fig. [28]) shows an unstable window of 
Pi £ (0.23, 0.5)w s . For 0.23u s < p\ < 0.5uj s , the S plot is above h(d) = 450000 and the converter is 
unstable with subharmonic oscillation. The unstable window of p\ is confirmed by time simulation. 
For p\ = 0.2u; s , the ripple size of y is small, and y°(t) is stable (Fig. [29]) . For p\ = 0.5w s , the 
ripple size of y is larger, and y°(t) is unstable (Fig. l27l) . For p\ = 0.6w s , the ripple size of y is even 
larger, but y°(t) is stable (Fig. [301) . Comparing Figs. [29][30l the ripple size of y°(t) is unrelated to 
subharmonic oscillation. This shows a counter-example for the hypothesis proposed in [U] that the 
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Figure 26: Curves of (| 13|) (solid line) and (|3ip (dashed line), the intersection is the subharmonic 
oscillation condition 




ripple size of y (t) is related to subharmonic oscillation. 

The unstable window of p\ is also confirmed by the sampled-data pole trajectories. The sampled- 
data pole trajectories for 0.1w s < p\ < 0.6oj s are shown in Fig. [3TJ Three poles are fixed around 
0.9485, 0.8853, and 0.51. A pole leaves the unit circle through -1 when p± = 0.23w s , and enters the 
unit circle when p± = 0.5iv s . This explains exactly the unstable window of p\. In this example, the 
unstable window of p\ is verified by three different approaches: time simulation, the S plot, and 
the sampled-data pole trajectories. □ 

Case II: p\ < u s /W. 

Similar to ACMC and based on simulations, for u p < uj s /10, the boundary condition (|16j) is still a 
good approximation. Here, CB\\ = and CA\Bn = — p\p2pR c K c / Z1Z2L = —pp2K c / k z , then the 
boundary condition (|16p becomes 

V* « WhKz ( I ) (54) 
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Figure 28: The intersection of the S plot (solid line) and h{d) (dashed line) shows the unstable 
window of p\ E (0.23, 0.5)o/ s . 




Figure 29: Stable T-periodic y°(t), pi = 0.2a/. 



6 Harmonic Balance Analysis of the Buck Converter 

A brief review of harmonic balance (HB) analysis is presented. The HB analysis is useful for 
both steady-state analysis and small-signal analysis of nonlinear systems [10^ [T8l [191 I20j . Since 
subharmonic oscillation of the boost converter under PVMC in continuous conduction mode does 
not occur as discussed in Sec. I4.3| without loss of generality, only the buck converter is considered. 

Consider a buck converter power stage, with a control-to-output (D-to-v ) transfer function 
G v d(s). In the converter, there is an ON switch and an OFF switch (sometimes substituted by 
a diode). Let the voltage across the OFF switch (or the diode) be Vd (as shown in Fig. [5j for 
example). The waveform of Vd(t) is a square wave with the high voltage at v s and the low voltage 
at 0, which can be represented by Fourier series (harmonics). 

In the converter, some parts are linear (from Vd to y) and some are nonlinear (from y to vj). 
Let the Vd-tov transfer function be G v (s). One has [JJ p. 470] 

g v (s) = 9}M = ^RcC + i (55) 
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Figure 30: Stable T-periodic y (t) with a larger ripple, p\ = 0.6a;, 




Figure 31: Sampled-data pole trajectories for 0.1u; s < pi < 0.6u; s . 



Let the compensator transfer function (from v D to —y (negative sign due to the negative feedback)) 
be G c (s). Let the transfer function from v& to — y be G(s). Then, for VMC, G(s) = G c (s)G v (s) = 
G c (s)G v d(s)/v s . Note that, for CMC, G(s) has a similar form as discussed later. Assume the PWM 
modulator has a gain 1/Vh- Let the loop gain be T(s), then T(s) = G c (s)G vc [(s)/Vh. The gain 
G(s) is proportional to the loop gain by 

G(s) = ^T(s) (56) 
v s 

The intersection of h(t) with the T-periodic solution y°(t) = Cx°(t) + Du determines the duty 
cycle and hence the waveform of Vd(t). By "balancing" the equation y°(t) = h(t) (written in 
Fourier series form) at the switching instants, conditions for existence of periodic solutions and 
subharmonic oscillation can be derived. Let Re denote taking the real part of a complex number. 
Based on [10} [20] , a necessary and sufficient condition for the occurrence of subharmonic oscillation 
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in a buck converter (with the trailing edge modulation) is v s < V* , where 



V 



2Re 



E [(1 " e^ D )G(jkco s ) - G(j(k - §) w ,)] 

k=l 



(57) 



expressed in terms of Vh, D, and the transfer function G(s) evaluated at half the switching frequency 
and its harmonics. Since both fjl3j) and (|57|) are exact conditions for the occurrence of subharmonic 
oscillation, it can be shown that f)13|) and (|57p are equivalent, but expressed in different forms. 
The boundary condition ([57]) also leads to a different expression of the S plot. Prom ([571) . 



1 

2^/ s Re[^((l - e^ kD )G(ju s k) - G(ju s (k - -)))] = h(d) 

k=l 



(58) 



where the left side of (j58|) is the S plot expressed in terms of harmonics. One can prove that (j58p 
is equivalent to (|12p . where the S plot is expressed in terms of matrices. 

Generally, G v (s), G c (s) and thus G(s) = G c (s)G v (s) are low-pass filters. The denominator of 
(|57p can be approximated by the term that involves G(s) with the smallest argument. Therefore, 
([57]) becomes 

Vh (59) 



V* 



2Re 



(1_ 6,^)0-^) _ G (i^) 



The HB analysis is applied to the various control schemes discussed above. For switching 
frequency much larger than l/V LC and 1/RC, the condition (|57p can be further simplified and 
expressed in closed forms. The derived closed-form boundary conditions are almost identical to 
those based on the sampled-data slope-based analysis discussed above, further corroborating the 
accuracy of the derived conditions. 



6.1 Proportional Voltage Mode Control (PVMC) 

Here, G c (s) = k p and G(s) = G c (s)G v (s) = k p G v (s). Based on the facts that for < D < 1 



then from ([55]) and ([5 

V* 



E 

k=l 



1 - cos(2irkD) 



k 2 



tt 2 L>(1 - D) 



E 



1^2 



E 

k=l 



sin(2vrfcL>) 
k 



TT' 



(60) 
(61) 
(62) 



V h LCuj 



k p p 
AV h LC 



r, r-v / j-, i sin(2-7rfcD) \ . 1 cos(2irkD)— 1 • 
1 1^ \-U s ti c ^{ g J + 77— TTa H p . 



fe=l 



fcpT 2 /) l iM(fl-i)|l-2i} + 2^ 



(63) 
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Figure 32: The intersection of the two curves (|63|) (solid line) and (|3ip (dashed line) shows V* = 24.5 



which is similar to f)28f) (based on the slope-based analysis) for R c <C -R (which results in p ~ 1) 
and R\ <C -L/C Since the harmonic balance analyzes the converter in the frequency domain, it 
has some advantages over the slope-based analysis. For example, it can show how the loading R 
affects the subharmonic oscillation as discussed later. 




Example 6. Consider a widely studied buck converter in [23]. The converter parameters are 
T = 400 fia, L = 20 mH, C = 47 fjF, R = 22 ft, V h = 4.4, and k p = 8.4. To transform from 
the leading-edge modulation to the generally used trailing-edge modulation and to remove the 
ramp offset in [24], let v r = 12.276, which can be proved to generate equivalent results as in [24]. 
Subharmonic oscillation is known to occur with V s * = 24.5. The intersection of the two curves (|63p 
and J3D also shows V* = 24.5 in Fig. [32] □ 



Dependence on loading resistance R. 

In ([63]) . V s * is almost independent of R because p = R/(R + R c ) « 1 for R c <C R- For practical 
converters, co s > 1/RC, and the frequency ratio r := l/RCu s is small. Then, (f6"3"|) is a good 
approximation. For large r, the dependence of V* on R becomes significant, which is studied next. 
Without loss of generality, assume that R c /R is small. 
From ([55]) and ([5j 



VhLClo 2 1^ 1 co S (27rkD)-l-lsm(27rkD) 

cos(27rD) - 1 - rsin(27rl?) 1 cos(4ttD) - 1 - § sin(47rD) 



i+r 2 1 + r 2 f+r 2 4 + r 



The next example illustrates the dependence on R. 



Example 7. [Dependence on R) Consider again Example [6] which has r = 0.062. Now decrease 
the loading resistance to R = 10 ft, then r = 0.136 becomes significant. Dependence on R is 
reported in [331 Fig- 7], which shows Vy = 26.8. The approximate curve (|65|) and the exact curve 
([57]) (also the same as ([T3]) ) are shown in Fig. [33] Based on ([57]) . = 26.8, agreed exactly with 
[33] . Based on ([65]) . V^* = 28, also agreed closely with [33]. □ 
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Figure 33: The intersection of the two curves (|57p (solid line) and (|3ip (dashed line) shows V* = 
26.8. The intersection of the two curves ([63]) (dotted line) and (f3Tj) (dashed line) shows V* = 28 



6.2 CMC with the Voltage Loop Open 

In CMC with the voltage loop open, the controlled output is II- Similar to ([55]) . the v^-to-iL 
transfer function is [TJ p. 470] 



a+A)gi+i 

LC(l + ^)s 2 + {^ + R c C)s + l 



Gi(s) := — , " * th , g— ■ - (66) 



Since no extra compensator (except the compensating ramp h(t)) is added in the current loop, 
G c (s) = 1, and one has G(s) := G c (s)Gi(s) = Gi(s) for CMC. From §3) and @2}, 

y* ~ -VhL = VhL 

S ~ T ^ sin( 27 rfcD) T s p(D - |) 
k=l 

For 72 c <C i?, ([671) is close to (f39l) (based on the slope-based analysis). 
6.3 CMC with the Voltage Loop Closed 

Here, there are two feedback loops. Let G c (s) be the voltage loop compensator. In the s-domain, 
V = i c ~ ib = G c (s)(v r - v ) -ih = G c (0)v r - (G c (s)G v (s) + Gi(s))v d . Then, 

G{s) = G c (s)G v (s) + G t (s) (68) 

Without loss of generality, let G c (s) = k p . Then, based on (|55p and (|66p . the boundary condition 
(|57p rearranged in terms of k p is 

m a L i T(_L I Rc\ nil 

fr * PS ti s + 4' RC ' L U + 2 (f , q) 

V T {1 ^2D 1) - + f) + (D - l)R c [ } 

which is similar to (|4ip (based on the slope-based analysis). 

Example 8. Consider again Example [2 Based on (|69p . one has fc* = 229, agreed closely with the 
simulation results (k* = 237). The small discrepancy between the HB analysis and the simulation 
results is due to the error of the duty cycle as discussed in Example El □ 
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6.4 Average Current Mode Control (ACMC) with a Type II Compensator 



The analysis of the ACMC buck converter in [12] is presented here for completeness. The de- 
nominator of (|57p can be approximated by the term that involves G with the smallest argument, 
then 

V* « — : (70) 

2Ke[G(ju s ) - G(^)] 

For an ACMC buck converter with oo s 3> 1/VLC and 1/RC, (ITOj) can be simplified as 

v; « Vj ^m (7i) 

where the frequency ratio 9 := lu p /uj s and 

m = ^±g£±jg (72) 

The function ip{6) = ?Jj(uj p /lj s ) has a minimum value of 5 at uj p = 0.38o; s . From ()T1 1) . it implies 
that larger values of ip, Vh, L, co z , and u a , or smaller values of R s and K c would lead to a larger 
stable operating range of the source voltage v s . This gives insight on how these parameters affect 
subharmonic oscillation. A larger value of K c /u z leads to larger mid-frequency gain of the current- 
loop compensator and larger crossover frequency of the current loop [341 135|. However, that also 
leads to instability. 

A design guideline proposed in suggests that the value of ui p is chosen between 0.33w s and 
0.5iv s . From Fig. [20] this choice of co p results in smaller V^* and hence smaller stable operating 
range of the source voltage. Setting uj p = u s instead, for example, will have larger operating range 
of the source voltage. From (fTT|) . V* is a function of u p . As u p affects stability nonlinearly, it is 
desirable to remove such effect (dependence on u p ). Let 

T/*| — T/*l ~ 5VhLu; z f s 

V s \min ■— V s lfl=0.38 ~ K 

which is not a function of u p . If v s < V*\ min , the converter is free from period-doubling for all u p . 
If v s > V*\ min , then there exists an unstable (period-doubling) window of uj p , and the converter is 
stable by choosing a larger value of oj p . Note that V*\ min is the minimum critical source voltage 
for all Up. It is possible that for a particular lo p and v s > V*\ min , the converter is still stable. But 
to eliminate the dependence on oj p , V*\ min is the maximum allowable source voltage with stability 
for all (jj p . 

Another stability guideline proposed in [28] suggests, expressed in terms of v s , 

2 1 V h Lu z f s W h Lu z f a , . 

v., < mm , — < (74 

h-D'D 1 R S K C R 8 K C V ' 



This stability guideline is more conservative than (|73p . The converter may be stretched as in ([73 
to achieve better dynamics. 

6.5 VMC with a Type III Compensator 

With the compensator (|53|) . which cancels some poles and zeros of the power stage, one has 

K 



G(s) = G c (s)G v (s) w ^ " 2sN (75) 
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Figure 34: Plot of (fi(D), solid line for (|78p and dashed line for more accurate (|77p 



which depends on uj s , and the compensator parameters k z and K c . From (|57p and ()75|) . the 
subharmonic oscillation is avoided if and only if 

v s < V* = ^^<t>(D) (76) 



where <p(D) is 

4>{D) 



Re 



' k(j-2k) (fc-i)(j-2fe+l) 



(77) 



A smaller k z or a larger K c would result in a smaller V* . Generally, uj s is large, and it may 
seem that the condition (|76p is generally met and the subharmonic oscillation is avoided. However, 
a large value of K c (in the order of u s ) may result in v s > V* and therefore subharmonic oscillation 
occurs. If (|59p . instead of (|57p . is used, the following approximate expression for (j){D) is obtained 

3 + 2 cos(z7rL>) — sm(27rL>) 
Both (J77[) and ([7S]) are shown in Fig. [Ml 

Example 9. {Effects of zero location) Consider again Example [4] where Z2 = l/2\/ LC (or 
k z = 1/2), and one has V* = 16. Now change to 1/VlC (or = 1) to see its effects. Simulation, 
confirmed with the exact sampled-data analysis, shows that the subharmonic oscillation occurs at 
V* = 23.9 V (and D ~ 0.138). Compared with Example [U one sees that a larger k z results in 
wider operating range of the source voltage. This effect can be clearly seen in (I76p which shows 
that a larger k z results in a larger V* . The value of V^* is not linearly proportional to k z because 
a different v s results in a different duty cycle to regulate to the same output voltage, and thus a 
different value of (fi(D). 

The subharmonic oscillation can be accurately predicted by the harmonic balance analysis (or 
equivalently the slope-based analysis). The plots of (|76j) (same as ([T3]) ) and (f3T1) are shown in 
Fig. [351 intersecting exactly at (D, V*) = (0.138, 23.9). Compared with Fig. [26l the curve ([76}, but 
not V*, is multiplied by two because k z is also multiplied by two. □ 
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Figure 35: Curves of (|76p (solid line) and (|3ip (dashed line), the intersection is the subharmonic 
oscillation condition 



7 Prediction of Subharmonic Oscillation Based on Loop Gain 



Since the gain G(s) is proportional to the loop gain T(s) by (|56p . and the subharmonic oscillation 
is avoided if and only if v s < V* , where V^* is given in (|57p . which directly (by simple algebra) 
leads to the following theorem. This theorem shows that a limitation on the loop gain is required 
to avoid the subharmonic oscillation. 



Theorem 1 Given a closed-loop buck converter with a loop gain T(s) 
avoided if and only if 



subharmonic oscillation is 



Re 



,fc=i 



J2k7TD 



)T(jkuj s 



T(j(k - 



1 

<2 



(79) 



The traditional loop gain analysis to design a converter with enough gain and phase margins 
only helps to avoid the Neimark-Sacker and the saddle- node bifurcations, but not the subharmonic 
oscillation, unless the system dimension is increased by considering the sampling effect as in [HE1E]. 
This theorem, putting an additional limitation on the familiar loop gain, helps to fill this gap. 

Using only one term in the series summation, (|79p becomes 



Re 



< 



1 



(80) 



As the following analysis indicates, a conservative design to make the converter free from the 
subharmonic oscillation for all duty cycles is by setting D = 1 in (|80p . Then (|80j) becomes 



Re 



T(if) 



> 



(81) 



which can be expressed in terms of Nyquist plot. In the Nyquist plot, the real part of the loop gain 
at half the switching frequency needs to be greater than —1/2 to avoid the subharmonic oscillation 
for all duty cycles. Therefore, the Nyquist plot is not only helpful to predict the Neimark-Sacker 
and the saddle-node bifurcations, it is also helpful to predict the subharmonic oscillation based on 
([HI]). Since the Nyquist plot is widely used, (f8"Tj) is helpful to design a stable converter. However, 
it should be noted that (|8ip is only an approximate condition, and the exact condition is (|79p . 
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"HB plot": a Nyquist-like plot in the complex plane. 

Note that the exact condition (f79|) is a function of D, uj s , and the loop gain T(s), where T(s) is 
further a function of v s , Vh, and the power stage and compensator parameters. Without loss of 
generality, the following design procedure for the line regulation design is proposed. In the line 
regulation design, one wants to find the stable operating range of v s . Given an desired output 
voltage, this is equivalent to find the stable operating range of D. Similar design procedures can 
be applied to determine the operating range of other parameters. Let 

oo 

H(D) := £(1 - e^ D )T{jku s ) - T(j(k - -)^ s ) (82) 

Then, the boundary condition (|79|) to avoid the subharmonic oscillation becomes 

Re [#(£>)] <± (83) 

For designation purpose, H{D) is called an HB plot because it is similar to the Nyquist plot to 
facilitate the converter design. Generally the feedback gain is large for a practical buck converter, 
then v s ~ v r /D. This value of v s may be used for the HB plot. Given a desired range of D, one 
can plot H(D) according to (|82p . Those values D such that the HB plot is on the left of 1/2 in the 
complex plane are the stable operating range of D to avoid the subharmonic oscillation. 



"M plot": a similar HB plot in matrix representation. 

Since both (113j) and (|57p are exact conditions for occurrence of subharmonic oscillation, they are 
actually identical. From (fTgj) and (E 



Re 



oo 1 

£[(1 _ e ^ D )G(jku s ) - G(j(k - -K)] 



,k=l 



|C[(I _ e AiTyi^ Al d _/) + (/ + e^y^Bu 



(84) 

Based on (|13p and (|28p . and similar to the HB plot in the complex plane, define a new "M plot" in 
the matrix form and as a function of D in the real domain, 

M{D) := ^C[(J- e AlT )- 1 (e AlTD -J) + (J + e AlT )- 1 ] J B 1 i (85) 

kr>T Da p ^:H,rC . , 1 . j C >, „Oxn j - 

* i^d [ -^ (D -2 )+p{1 -^L- ){1 - 2D+2D )] (86) 



where v s ~ v r /D for a large feedback gain (which is true for most practical converters). Note 
that (|85p is an exact representation, and (|86p is an approximate representation, which is generally 
accurate if the controller poles are smaller than w s /10 as discussed previously. Then, from Theorem 
1, (ESI) and (1841), one has 



Corollary 1 Given a closed-loop buck converter as in Fig. [H subharmonic oscillation is avoided if 
and only if 

M(D) < 1 (87) 



Actually, this corollary can be also proved easily from (I13p (or (|28p ). by dividing both sides of (|13p 
(or (|28p ) by the right side (to normalize the right side to one). Both (|85j) and (j 13j) are equivalent 
exact boundary conditions but expressed in different forms. 
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Figure 36: HB plot in the complex plane showing the stable operating range of D = [0.34, 0.89]. 
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Figure 37: M plot M(D) showing the stable operating range of D = [0.34,0.89]. 



Example 10. {Line regulation: determine the stable operating ranges of D and v s , based on the 
HB plot and the M plot.) Consider again Example [U but with R c = 2 mfi. In Example Q3 it is 
shown that the stable operating range of D is [0.34,0.89]. Here, using the HB plot and the M plot 
will confirm the same result. 

An HB plot for D = [0.25,0.99] is shown in Fig. [36l intersecting with the vertical line at 1/2 
in the complex plane at D = 0.34 and 0.89. An M plot for D = [0.25, 0.99] is shown in Fig. E3 
intersecting with the horizontal line at 1 at D = 0.34 and 0.89. Both the plots agrees exactly 
with the analysis in Example [TJ Given an output voltage around v r = 4, the corresponding stable 
operating range of v s free from the subharmonic oscillation is [4/0.89,4/0.34] = [4.48, 11.85]. □ 

Example 11. (Application of Theorem 1.) Consider again Example [61 but with R = 2 Q, and 
v s = 50 V. Using the exact sampled-data analysis, D = 0.243, x°(0) = (5.9867,12.0753)', and 
x °(d) = (6.1711,12.1486)'. The eigenvalues of $ are -0.0336 and -0.4222, and they are stable. 
The steady-state simulation is shown in Fig. [38l also indicating that the converter is stable. From 
flggD , #(0.243) = 0.1390 + 0.8867j and Re [if (0.243)] < 1/2. Based on Theorem 1, the converter 
is stable, agreed with the sampled-data analysis and the simulation. From (|13p. subharmonic 
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Figure 38: Plots of y(t) and 2l(£) (solid lines), and h(t) (dashed line), R = 2 and v s = 50 



oscillation occurs when v s > 82.9. 



□ 



8 Subharmonic Oscillation Caused by a Large Crossover Frequency 

A rule of thumb to select the crossover frequency uj c is by setting w c < u s /5 [22 p. 250]. No 
theoretical explanation is given. The following analysis will give a theoretical explanation and puts 
a limit on the crossover frequency to avoid the subharmonic oscillation. Since the compensator 
poles and zeros are selected based on the power stage parameters, such as L, C, and R c , the 
selected poles and zeros may result in a large crossover frequency, and the subharmonic oscillation 
occurs. 

The loop gain analysis in Theorem 1 can lead to a condition based on the crossover frequency. 
Without loss of generality, only two control schemes are considered. Similar analysis can be applied 
to other control schemes. 

8.1 Buck Converter with a Type III Compensator 

Case I: p\ = oo s /2 > w s /10. 

First, a closed-form crossover frequency is derived. The loop gain (after some poles and zeros being 
canceled) for the frequency below oj s /2 is approximately 



From (|88|) . the crossover frequency lo c is around K c v s / n z Vh- In Example it is shown that a 
smaller k z leads to a smaller operating range of v s . From (|88p . a smaller k z leads to a larger lo c , 
and the subharmonic oscillation is more likely to occur as discussed next. With the estimated 
crossover frequency u^ st := K c v s / K z Vh, rearranging (|76p leads to the following theorem. 

Theorem 2 In a buck converter with a type III compensator A53\) , the subharmonic oscillation is 
avoided if and only if 







where 4>(D), shown in (77), varies from 0.694 t° (Fi>9- \^4\) - 
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A larger K c v s or a smaller n z Vh leads to a larger lo c to a point such that (|89p is violated and the 
subharmonic oscillation occurs. Since the minimum of <f)(D) is 0.694, a rule of thumb to avoid 
the subharmonic oscillation for all duty cycles is lo c < 0.347w s . This theorem gives a theoretical 
explanation on the traditional wisdom on setting oo c < u s /5. In Example [U with k z = 0.5, the 
subharmonic oscillation occurs when uj c = 1.08 x 10 6 = 0.58cj s (shown in Fig. |25|) . where uj c is too 
large. 

From Fig. [341 the maximum of (J)(D) occurs around D = 0.4 with a value of 2.89. Therefore, 
around D = 0.4, only u c < 1.44o; s is required. The limitation on the crossover frequency for 
operating at this duty cycle may be relaxed. However, to have a wider operating range of D, 
oj c < 0.347o; s is probably a better design (although it may be conservative in some situations). 

Case II: pi < co s /W. 

As the analysis above, the approximate crossover frequency is derived first. The loop gain for 
frequency below oj s /2 is approximately 

Solving \T(jw c )\ = 1 gives 



w ^~Y + p + ^n (91) 

For p2 < w s /10, the approximate condition ([54"]) is close to the exact condition (fT3|) as discussed 
above. The curve of the condition (j54[) has a minimum at D = 1 (or D = 0) with a value of 



T 2 p 2 K c 



V:\ min = V:\ D=l = ^r (92) 



A conservative design to avoid subharmonic oscillation for all duty cycles is v s < V*\ a . Then, 
(I9TD leads to 



Pi , JP2 , w s 



For £>2 = w s /10, (|93p leads to uj c < 0.27u s , which is required to avoid the subharmonic oscillation 
for all duty cycles. For p 2 <C oj s , based on ([93]) . oj c /^s < 1/vr is required to avoid the subharmonic 
oscillation for all duty cycles. 

8.2 ACMC with a Type II Compensator 

First, the approximate crossover frequency is derived. Here, the loop gain at high frequency is 
rj,, v v s R S r , f . n i v v s Rs %7 Cs v s R s K c 

m = — GcWGlW - = —j. _ _ (94) 

Solving \T(jw c )\ = 1 gives 



w P , u l , uj p v s K c R s 

w ^-Y + H + ^l- (95) 
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Without loss of generality, only the case for u p < cj s /10 is considered. For cj p < uj s /W, the 
approximate condition (|46p is close to the exact condition (|13p as shown in Fig. [2TJ The curve of 
the condition (|46p has a minimum at D = 1 (or D = 0) with a value of 



1 *a 



= v* 



4V h oj z L 



(96) 



s \rnin 



s \D=1 ~ 



T 2 K c R s u p 



A conservative design to avoid subharmonic oscillation for all duty cycles is v s < V*\ min . Then, 
(195]) leads to 



For ui p = lu s /10, based on ()97|) . uj c < 0.27o; s is required to avoid subharmonic oscillation for all duty 
cycles. For uj p <C oj s , based on ([93j) . uj c /uj s < 1 /it is required to avoid the subharmonic oscillation 
for all duty cycles. The analysis is almost identical as for the type III compensator discussed above. 

9 Conclusion and Contributions 

Based on a subharmonic oscillation boundary condition published in a PhD thesis (but not else- 
where) more than a decade ago, extended design-oriented boundary conditions are derived for 
general switching DC-DC converters expressed by a unified VMC/CMC block diagram model. Un- 
der the unified framework, the boundary conditions for both VMC and CMC have similar forms. 
One does not need to analyze VMC and CMC separately. The boundary conditions are expressed 
in terms of signal slopes. The well known slope-based instability criterion for CMC becomes a 
special case in this unified modeling approach. The derived sloped-based boundary conditions are 
exact, and can be further simplified in various approximate closed forms to facilitate the converter 
design. The sloped-based boundary conditions are applied to analyze various VMC/CMC control 
schemes of varying complexity. Among the various control schemes, those popular type II or type 
III controllers (with integrators included) are also analyzed. 

Harmonic balance analysis is also applied to derive the exact subharmonic oscillation boundary 
conditions, which are equivalent to the sloped-based boundary conditions. It is also applied to 
analyze five different VMC/CMC control schemes. Based on the harmonic balance analysis, a new 
"HB plot" in the complex plane, similar to the Nyquist plot, is proposed to accurately predict the 
occurrence of the subharmonic oscillation. Another equivalent "M plot" in the real domain is also 
proposed, which also accurately predicts the occurrence of the subharmonic oscillation. Based on 
the harmonic balance analysis, the relation between the crossover frequency and the subharmonic 
oscillation can be clearly seen. 

Both the sloped-based analysis and the harmonic balance analysis complement each other. 
The slope-based analysis analyzes the converter in the time domain and expresses the boundary 
conditions in matrix forms. The harmonic balance analysis expresses the boundary conditions 
in terms of (sub)harmonics of the switching frequency. Since most power stage and compensator 
transfer functions are low-pass filters, the harmonic balance analysis is particularly useful to analyze 
the converter with the controller poles close to the switching frequency. 

Both the slope-based analysis and the harmonic balance analysis complement the traditional 
average analysis. The average analysis can predict the Niemark-S acker bifurcation and the saddle- 
node bifurcation, but generally not the subharmonic oscillation, unless the sampling effects are 
considered. In Example [21 the converter has a phase margin of 129 degrees based on an average 
model, but the subharmonic oscillations still occur. The subharmonic oscillation in all examples 
presented can be accurately predicted by the slope-based analysis or the harmonic balance analysis. 




(97) 
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The contributions of this paper are as follows. First, the following sixteen boundary conditions 
are reported here for the first time to the author's knowledge. Some of the conditions are equivalent 
to each other. These derived boundary conditions greatly facilitate future converter design to avoid 
subharmonic oscillation. Since there are many equations reported in this paper, the key equations 
are identified here for easy reference. The most general boundary condition for the subharmonic 
oscillation is Q. It is an exact condition expressed in terms of signal slopes. For the boost converter, 
© becomes (|18p . which is also an exact condition. For the buck converter, © becomes (|12p . which 
is equivalent to various forms: (|13[) in terms of v s , ()22[) in terms of the S plot, (j79|) in terms of the 
loop gain, (|83p in terms of the HB plot, and (jS7J) in terms of the M plot. These boundary conditions 
can be further simplified to various approximate forms: (|15p in terms of signal slopes, (|16p in terms 
of v s , and (j80|) in terms of the loop gain. For those control schemes presented, the boundary 
conditions can be further simplified to more compact forms, such as (|28p for PVMC and V 2 control 
(with ESR considered), (f37"|) for CMC with the voltage loop open (with ESR considered), (|4"l~]l for 
CMC with the voltage loop closed, (|46|) for ACMC with a type II controller and small controller 
poles, and (|54p for VMC with a type III controller and small controller poles. 

Second, all the boundary conditions are in closed forms, one can see the effects of various con- 
verter parameters on the stability. In some control schemes, the effects of ESR, loading resistance, 
feedback gain, and pole/zero locations are further analyzed. For example, in CMC, ESR extends 
the stability beyond D = 1/2 as shown in (|38p . Also, for CMC or VMC, ESR causes the boundary 
conditions to have both terms D — 1/2 and 1 — 2D + 2D 2 . 

Third, although only some control schemes are analyzed, this paper also provides a methodology 
based on the general boundary condition to analyze other converters with different control schemes 
to predict the subharmonic oscillation. As discussed in the Introduction, for a converter under a 
particular control scheme, once the converter is expressed in terms of the unified VMC/CMC block 
diagram model shown in Fig. 1, the boundary condition for that particular control scheme can be 
readily obtained. 

Fourth, under this general methodology, those recently reported boundary conditions for PI 
control [9] and V 2 control [UJ become special cases in view of the general boundary condition. In 
consideration of recent development to derive more boundary conditions for other converters or 
control schemes, the boundary conditions derived here greatly facilitate future research. 

Fifth, new Nyquist-like design-oriented plots (such as the S, HB, M plots) to predict the occur- 
rence of subharmonic oscillation are proposed. The S plot also shows the required ramp slope to 
prevent the occurrence of subharmonic oscillation. 

Sixth, closed- form conditions (such as (|59"1) . (J9"3"j) and (f9"T|) ) relating the crossover frequency to 
the subharmonic oscillation are derived. 
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